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Abstract. There are many reasons to believe the present mass 
density of the universe is dominated by a weakly interacting 
massive particle (wimp), a fossil relic of the early universe. The- 
oretical ideas and experimental efforts have focused mostly on 
production and detection of thermal relics, with mass typically 
in the range a few GeV to a hundred GeV. Here, I will review 
scenarios for production of nonthermal dark matter. Since the 
masses of the nonthermal wimps are in the range 10 12 to 10 16 
GeV, much larger than the mass of thermal wimpy WIMPS, they 
may be referred to as wimpzillas. In searches for dark matter 
it may be well to remember that "size does matter." 
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1. Introduction 



There is conclusive evidence that the dominant component of the matter 
density in the universe is dark. The most striking indication of the ex- 
istence of dark matter is the dynamical motions of astronomical objects. 
Observations of flat rotation curves for spiral galaxies Q indicates that the 
dark component of galactic halos is about ten times the luminous com- 
ponent. Dynamical evidence for DM in galaxy clusters from the velocity 
dispersion of individual galaxies, as well as from the large x-ray tempera- 
tures of clusters, is also compelling Bulk flows, as well as the peculiar 
motion of our own local group, also implies a universe dominated by dark 
matter ||. 

The mass of galaxy clusters inferred by their gravitational lensing of 
background images is consistent with the large dark-to-visible mass ratios 
determined by dynamical methods 

There is also compelling evidence that the bulk of the dark component 
must be nonbaryonic. The present baryonic density is restricted by big- 
bang nucleosynthesis to be less than that inferred by the methods discussed 
above || . The theory of structure formation from the gravitational instabil- 
ity of small initial seed inhomogeneities requires a significant nonbaryonic 
component to the mass density 

In terms of the critical density, pc = SHqM^/Stt = 1.88 x 10~ 29 g cm" 3 
with Hubble constant Hq = lOO/i kmsec _1 Mpc _1 and Planck mass Mp\, 
the dark-matter density inferred from dynamics is Odm = Pbm/pc ~ 0.3. 
In addition, the most natural inflation models predict a flat universe, i.e., 
Qq = 1, while standard big-bang nucleosynthesis implies that ordinary 
baryonic matter can contribute at most 10% to f2 - This means that about 
90% of the matter in our universe may be dark. 



2. Thermal Relics — Wimpy WIMPS 

It is usually assumed that the dark matter consists of a species of a new, yet 
undiscovered, massive particle, traditionally denoted by X. It is also often 
assumed that the dark matter is a thermal relic, i.e., it was in chemical 
equilibrium in the early universe. 

A thermal relic is assumed to be in local thermodynamic equilibrium 
(lte) at early times. The equilibrium abundance of a particle, say relative 
to the entropy density, depends upon the ratio of the mass of the particle to 
the temperature. Define the variable Y = nx/s, where nx is the number 
density of WIMP X with mass Mx, and s ~ T 3 is the entropy density. 
The equilibrium value of Y , Yeq, is proportional to exp(— x) for x ^> 1, 
while Yeq ~ constant for x <C 1, where x = Mx/T. 

2 








Y 



«-Y, 



EQ 



-20 



1 



10' 



10 5 " 



10 3 



x = M x /T 



Figure 1. A thermal relic starts in lte at T > Mx- When the rates 
keeping the relic in chemical equilibrium become smaller than the expansion 
rate, the density of the relic relative to the entropy density freezes out. 

A particle will track its equilibrium abundance as long as reactions 
which keep the particle in chemical equilibrium can proceed rapidly enough. 
Here, rapidly enough means on a timescale more rapid than the expansion 
rate of the universe, H . When the reaction rate becomes smaller than the 
expansion rate, then the particle can no longer track its equilibrium value, 
and thereafter Y is constant. When this occurs the particle is said to be 
"frozen out." A schematic illustration of this is given in Fig. [I]. 

The more strongly interacting the particle, the longer it stays in LTE, 
and the smaller its eventual freeze-out abundance. Conversely, the more 
weakly interacting the particle, the larger its present abundance. The 
freeze-out value of Y is related to the mass of the particle and its anni- 
hilation cross section (here characterized by ctq) by |tJ 



Since the contribution to is proportional to Mxnx, which in turn is 
proportional to MxY, the present contribution to fl from a thermal relic 
roughly is independent of its massj^] and depends only upon the annihila- 
tion cross section. The cross section that results in Q,xh 2 ~ 1 is of order 

1 To first approximation the relic dependence depends upon the mass only indirectly 
through the dependence of the annihilation cross section on the mass. 
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10 _37 cm 2 , of the order the weak scale. This is one of the attractions of 
thermal relics. The scale of the annihilation cross section is related to a 
known mass scale. 

The simple assumption that dark matter is a thermal relic is surpris- 
ingly restrictive. The largest the annihilation cross section can be is roughly 
M^ 2 . This implies that large-mass WIMPS would have such a small annihi- 
lation cross section that their present abundance would be too large. Thus 
one expects a maximum mass for a thermal WIMP, which turns out to be 
a few hundred TeV || . 

The standard lore is that the hunt for dark matter should concentrate on 
particles with mass of the order of the weak scale and with interaction with 
ordinary matter on the scale of the weak force. This has been the driving 
force behind the vast effort in dark matter direct detection described in 
this meeting by Cabrera jfj , Liubarsky jio) , Bernabei [[ll] , Ramachers O] , 
and Baudis |l3|j . 

In view of the unitarity argument, in order to consider thermal wim- 
PZILLAS, one must invoke, for example, late-time entropy production to 
dilute the abundance of these supermassive particles |l4|], rendering the 
scenario unattractive. 



3. Nonthermal Relics — WIMPZILLAS 

There are two necessary conditions for the wimpzilla scenario. First, the 
wimpzilla must be stable, or at least have a lifetime much greater than 
the age of the universe. This may result from, for instance, supersymmctric 
theories where the breaking of supersymmetry is communicated to ordinary 
sparticles via the usual gauge forces In particular, the secluded and 
the messenger sectors often have accidental symmetries analogous to baryon 
number. This means that the lightest particle in those sectors might be 
stable and very massive if supersymmetry is broken at a large scale JD|. 
Other natural candidates arise in theories with discrete gauge symmetries 
jl7j and in string theory and M theory Jl^, |l^] . 

It is useful here to note that wimpzilla decay might be able to account 
for ultra-high energy cosmic rays above the Greisen-Zatzepin-Kuzmin cut- 
off |2f|. A wimpy little thermal relic would be too light to do the job, 
a wimpzilla is needed. 

The second condition for a wimpzilla is that it must not have been in 
equilibrium when it froze out (i.e., it is not a thermal relic), otherwise fix h 2 
would be much larger than one. A sufficient condition for noncquilibrium is 
that the annihilation rate (per particle) must be smaller than the expansion 
rate: < H, where a\v\ is the annihilation rate times the M0ller flux 

factor, and H is the expansion rate. Conversely, if the dark matter was 
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created at some temperature T» and Clxh 2 < 1, then it is easy to show 
that it could not have attained equilibrium. To see this, assume X's were 
created in a radiation-dominated universe at temperature T». Then flxh 2 
is given by 

n x h 2 = n 7 /i 2 (T*/T )mxnx(T*)/ / 9 7 (T ! ,) , (2) 

where To is the present temperature. Using the fact that p 7 (T*) = 
H{T*)MpiT 2 , n x {T«)/H(T«) = (tt x /Cl 7 )T M P iT*/M x . One may safely 
take the limit a\v\ < M x , so n x (T*)a\v\/ 'H{T*) must be less than 
(£l x /Vtr^ToMpiT* I 'M x . Thus, the requirement for nonequilibrium is 



200TcV\ ( T, 



M x ) \M X 



< 1 . (3) 



This implies that if a nonrelativistic particle with M x > 200 TeV was cre- 
ated at T* < M x with a density low enough to result in Cl x ~ 1, then 
its abundance must have been so small that it never attained equilibrium. 
Therefore, if there is some way to create wimpzillas in the correct abun- 
dance to give fix ~ 1, nonequilibrium is automatic. 

Any wimpzilla production scenario must meet these two criteria. Be- 
fore turning to several wimpzilla production scenarios, it is useful to esti- 
mate the fraction of the total energy density of the universe in WIMPZILLAS 
at the time of their production that will eventually result in CI ~ 1 today. 

The most likely time for wimpzilla production is just after inflation. 
The first step in estimating the fraction of the energy density in wimpzil- 
las is to estimate the total energy density when the universe is "reheated" 
after inflation. 

Consider the calculation of the reheat temperature, denoted as Trh- 
The reheat temperature is calculated by assuming an instantaneous con- 
version of the energy density in the inflaton field into radiation when the 
decay width of the inflaton energy, T^, is equal to H 7 the expansion rate 
of the universe. 

The reheat temperature is calculated quite easily Q . After inflation the 
inflaton field executes coherent oscillations about the minimum of the po- 
tential. Averaged over several oscillations, the coherent oscillation energy 
density redshifts as matter: p^ cx a~ 3 , where a is the Robertson- Walker 
scale factor. If pi and ai denotes the total inflaton energy density and 
the scale factor at the initiation of coherent oscillations, then the Hubble 
expansion rate as a function of a is 



3 Mp t V a 



Equating H(a) and leads to an expression for ai/a. Now if all available 
coherent energy density is instantaneously converted into radiation at this 
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value of aj /a, one can define the reheat temperature by setting the coherent 
energy density, p$ = pi{ai/a) 3 7 equal to the radiation energy density, 
PR — ( 7I " 2 /30)5*T^ H , where g* is the effective number of relativistic degrees 
of freedom at temperature Trh- The result is 



(- 



90 \ 1/4 _/200^ 1/4 



(6) 



^-{S^-J Vm^i- (5) 

The limit from gravitino overproduction is Trh ~ 10 9 to 10 10 GcV. 

Now consider the WIMPZILLA density at reheating. Suppose the wim- 
PZILLA never attained lte and was nonrelativistic at the time of produc- 
tion. The usual quantity Vtxh 2 associated with the dark matter density 
today can be related to the dark matter density when it was produced. 
First write 

pxjtp) _ PxjtRn) I T RH \ 

Ar(*o) PR^Rh) V T J ' 

where pr denotes the energy density in radiation, px denotes the energy 
density in the dark matter, Trh is the reheat temperature, To is the temper- 
ature today, to denotes the time today, and tRH denotes the approximate 
time of reheating.^ To obtain px {tRH ) / Pr {tRH ) , one must determine when 
X particles are produced with respect to the completion of reheating and 
the effective equation of state between X production and the completion 
of reheating. 

At the end of inflation the universe may have a brief period of matter 
domination resulting either from the coherent oscillations phase of the in- 
flaton condensate or from the preheating phase |22| . If the X particles are 
produced at time t = t e when the de Sitter phase ends and the coherent 
oscillation period just begins, then both the X particle energy density and 
the inflaton energy density will redshift at approximately the same rate 
until reheating is completed and radiation domination begins. Hence, the 
ratio of energy densities preserved in this way until the time of radiation 
domination is 

PxjtRH) 8?r px{t e ) , s 

PR{t RH ) ~ 3 M^H^teY 1 > 

where Mpi w 10 19 GeV is the Planck mass and most of the energy density 
in the universe just before time tRH is presumed to turn into radiation. 
Thus, using Eq. |^, one may obtain an expression for the quantity Qx = 
Px {to)/ pc {to), where pc{h) = SH^M^/Sk and H = 100 h km sec" 1 
Mpc" 1 : 

2 More specifically, this is approximately the time at which the universe becomes 
radiation dominated after inflation. 
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Here 51 rIi 2 ss 4.31 x 10 -5 is the fraction of critical energy density in ra- 
diation today and nx is the density of X particles at the time when they 
were produced. 

Note that because the reheating temperature must be much greater 
than the temperature today (Trh/Tq > 4.2 x 10 14 ), in order to satisfy 
the cosmological bound flxh 2 < 1, the fraction of total wimpzilla energy 
density at the time when they were produced must be extremely small. 
One sees from Eq. | that tt x h 2 ~ 10 17 ' {T RH /10 9 GeV)(p x {t e )/ >(* e ))- It 
is indeed a very small fraction of the total energy density extracted in 

WIMPZILLAS. 

This means that if the wimpzilla is extremely massive, the challenge 
lies in creating very few of them. Gravitational production discussed in 
Section LI naturally gives the needed suppression. Note that if reheating 
occurs abruptly at the end of inflation, then the matter domination phase 
may be negligibly short and the radiation domination phase may follow 
immediately after the end of inflation. However, this does not change Eq. 



4. WIMPZILLA PRODUCTION 

4-1. Gravitational Production 

First consider the possibility that wimpzillas are produced in the tran- 
sition between an inflationary and a matter-dominated (or radiation- 
dominated) universe due to the "nonadiabatic" expansion of the back- 
ground spacetime acting on the vacuum quantum fluctuations p3[ . 

The distinguishing feature of this mechanism is the capability of gener- 
ating particles with mass of the order of the inflaton mass (usually much 
larger than the reheating temperature) even when the particles only in- 
teract extremely weakly (or not at all) with other particles and do not 
couple to the inflaton. They may still be produced in sufficient abun- 
dance to achieve critical density today due to the classical gravitational 
effect on the vacuum state at the end of inflation. More specifically, if 
0.04 < Mx/Hi < 2, where Hi ~ m$ ~ 10 13 GeV is the Hubble constant at 
the end of inflation {m<f, is the mass of the inflaton) , wimpzillas produced 
gravitationally can have a density today of the order of the critical density. 
This result is quite robust with respect to the "fine" details of the transi- 
tion between the inflationary phase and the matter-dominated phase, and 
independent of the coupling of the wimpzilla to any other particle. 

Conceptually, gravitational wimpzilla production is similar to the in- 
flationary generation of gravitational perturbations that seed the formation 
of large scale structures. In the usual scenarios, however, the quantum gen- 
eration of energy density fluctuations from inflation is associated with the 
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inflaton field that dominated the mass density of the universe, and not 
a generic, sub-dominant scalar field. Another difference is that the usual 
density fluctuations become larger than the Hubble radius, while most of 
the WIMPZILLA perturbations remain smaller than the Hubble radius. 

There are various inequivalent ways of calculating the particle produc- 
tion due to interaction of a classical gravitational field with the vacuum 
(see for example |24j ], [ p5[ , and (2(|). Here, I use the method of finding the 
Bogoliubov coefficient for the transformation between positive frequency 
modes defined at two different times. For Mx/Hj < 1 the results are quite 
insensitive to the differentiability or the fine details of the time dependence 
of the scale factor. For 0.04 < Mx/Hj < 2, all the dark matter needed for 
closure of the universe can be made gravitationally, quite independently of 
the details of the transition between the inflationary phase and the matter 
dominated phase. 

Start with the canonical quantization of the X field in an action of the 
form (with metric ds 2 = dt 2 — a 2 (t)dx 2 = a 2 (if) [drj 2 — dx 2 ] where 77 is 
conformal time) 

S = { dt { d 3 x^(x 2 - - M\X 2 - iRX 2 \ (9) 



2 

where R is the Ricci scalar. After transforming to conformal time coordi- 
nate, use the mode expansion 

/d 3 k r 1 
(27r) 3 /2a(??) (^ k - x + 4^K ik - x ] , (10) 

where because the creation and annihilation operators obey the commu- 
tator [afcj^aj^] = #( 3 )(ki — k 2 ), the h^s obey a normalization condition 
hkhfr — h! k h* k = i to satisfy the canonical field commutators (henceforth, 
all primes on functions of 77 refer to derivatives with respect to 77). The 
resulting mode equation is 

h'l(r,)+wl(r 1 )h k ( V )=0, (11) 

where 

w 2 k = k 2 + M\a 2 + (6£ - l)a"/a . (12) 

The parameter £ is 1/6 for conformal coupling and for minimal coupling. 
From now on, £ = 1/6 for simplicity but without much loss of generality. 
By a change in variable 77 — > k/a, one can rewrite the differential equation 
such that it depends only on H(rj), H'(r])/k, k/a(rf), and Mx- Hence, the 
parameters Hj and a/ correspond to the Hubble parameter and the scale 
factor evaluated at an arbitrary conformal time 77/, which can be taken to 
be the approximate time at which As are produced (i.e., r/i is the conformal 
time at the end of inflation). 
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Figure 2. The contribution of gravitationally produced wimpzillas to 
Ox/i 2 as a function of Mx/Hj. The shaded area is where thermalization 
may occur if the annihilation cross section is its maximum value. Also 
shown is the contribution assuming that the wimpzilla is present at the 
end of inflation with a temperature T = Hj/2tt. 




Figure 3. The evolution of the Bogoliubov coefficient with conformal time 
for several wavenumbers. r\ = r\i corresponds to the end of the inflationary 
era. 
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One may then rewrite Eq. |ll| as 

^) + (fc 2 + §|a 2 )^(?7)=0, (13) 

where fj — r/aiHi, a = a/ai, and k — k/(aiHi). For simplicity of nota- 
tion, drop all the tildes. This differential equation can be solved once the 
boundary conditions are supplied. 

The number density of the wimpzillas is found by a Bogoliubov trans- 
formation from the vacuum mode solution with the boundary condition at 
77 = 770 (the initial time at which the vacuum of the universe is determined) 
into the one with the boundary condition at 77 = 771 (any later time at 
which the particles are no longer being created). 770 will be taken to be 
—00 while 771 will be taken to be at +00. Defining the Bogoliubov transfor- 
mation as fi^iv) = a fc' l fc ( r 7) + Pkh* k Vo (r]) (the superscripts denote where 
the boundary condition is set), the energy density of produced particles is 

/ 1 \ 3 r°° dk - 
px{m) =M x n x ( m ) = M x Hll^—) j ^\(3 k \\ (14) 

where one should note that the number operator is defined at 771 while the 
quantum state (approximated to be the vacuum state) defined at 770 does 
not change in time in the Heisenberg representation. 

As one can see from Eq. [l|, the input parameter is Mx/Hi. One must 
also specify the behavior of 0(77) near the end of inflation. In Fig. |^ (from 
[ p3| ) , I show the resulting values of Vtxh 2 as a function of Mj /Hj assuming 
the evolution of the scale factor smoothly interpolates between exponen- 
tial expansion during inflation and either a matter-dominated universe or 
radiation-dominated universe. The peak at Mx/Hi ~ 1 is similar to the 
case presented in Ref. [|7|. As expected, for large Mx/Hi, the number 
density falls off faster than any inverse power of Mx/ 'Hi. 

Now most of the action occurs around the transition from inflation to 
the matter-dominated or radiation-dominated universe. This is shown in 
Fig. |. Also from Fig. | one can see that most of the particles are created 
with wavenumber of order Hi. 

To conclude, there is a significant mass range (O.li?/ to Hi, where 
Hi ~ 10 13 GeV) for which wimpzillas will have critical density today 
regardless of the fine details of the transition out of inflation. Because this 
production mechanism is inherent in the dynamics between the classical 
gravitational field and a quantum field, it needs no fine tuning of field 
couplings or any coupling to the inflaton field. However, only if the particles 
are stable (or sufficiently long lived) will these particles give contribution 
of the order of critical density. 
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4-2. Production during Reheating 

Another attractive origin for wimpzillas is during the defrosting phase 
after inflation. It is important to recall that it is not necessary to convert 
a significant fraction of the available energy into massive particles; in fact, 
it must be an infinitesimal amount. I now will discuss how particles of 
mass much greater than T R h may be created in the correct amount after 
inflation in reheating |2S|| . 

In one extreme is the assumption that the vacuum energy of inflation is 
immediately converted to radiation resulting in a reheat temperature T R h ■ 
In this case fix can be calculated by integrating the Boltzmann equation 
with initial condition N x = at T — T RH . One expects the X density to be 
suppressed by exp(— 2M X /T RH ); indeed, one finds fix ~ 1 for M X /T RH ~ 
25 + 0.51n(TO x (cr|u|)), in agreement with previous estimates pO that for 
Trh ~ 10 9 GeV, the wimpzilla mass would be about 2.5 x lO^GeV. 

A second (and more plausible) scenario is that reheating is not instanta- 
neous, but is the result of the slow decay of the inflaton field. The simplest 
way to envision this process is if the comoving energy density in the zero 
mode of the inflaton decays into normal particles, which then scatter and 
thermalizc to form a thermal background. It is usually assumed that the 
decay width of this process is the same as the decay width of a free inflaton 
field. 

There are two reasons to suspect that the inflaton decay width might 
be small. The requisite flatness of the inflaton potential suggests a weak 
coupling of the inflaton field to other fields since the potential is renormal- 
ized by the inflaton coupling to other fields |2S|] . However, this restriction 
may be evaded in supersymmetric theories where the nonrenormalization 
theorem ensures a cancelation between fields and their superpartners. A 
second reason to suspect weak coupling is that in local supersymmetric 
theories gravitinos are produced during reheating. Unless reheating is de- 
layed, gravitinos will be overproduced, leading to a large undesired entropy 
production when they decay after big-bang nucleosynthesis [ ]30| . 

It is simple to calculate the wimpzilla abundance in the slow reheating 
scenario. It will be important to keep in mind that what is commonly called 
the reheat temperature, T R h, is not the maximum temperature obtained 
after inflation. The maximum temperature is, in fact, much larger than 
T R h ■ The reheat temperature is best regarded as the temperature below 
which the universe expands as a radiation-dominated universe, with the 
scale factor decreasing as 1 ' 3 T~ 1 . In this regard it has a limited meaning 
[0, |3lj. One implication of this is that it is incorrect to assume that the 
maximum abundance of a massive particle species produced after inflation 
is suppressed by a factor of exp(—M/T RH ). 
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To estimate wimpzilla production in reheating, consider a model uni- 
verse with three components: inflaton field energy, p^, radiation energy 
density, pn, and wimpzilla energy density, px- Assume that the decay 
rate of the inflaton field energy density is T^. Also assume the wimpzilla 
lifetime is longer than any timescale in the problem (in fact it must be 
longer than the present age of the universe) . Finally, assume that the light 
degrees of freedom are in local thermodynamic equilibrium. 

With the above assumptions, the Boltzmann equations describing the 
redshift and interchange in the energy density among the different compo- 
nents is 



PR 



iHp^ - 
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px + 'iHpx 



m x 



m x 
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Px 



Px 



P E x Q 



(p E x Q 
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(15) 



where dot denotes time derivative. As already mentioned, (cr|v|) is the ther- 
mal average of the X annihilation cross section times the M0ller flux factor. 
The equilibrium energy density for the X particles, p x ®, is determined by 
the radiation temperature, T — (30pn/ir 2 g r ) 1 ^ 4 . 

It is useful to introduce two dimensionless constants, and ax, defined 
in terms of and (a\v\) as 

rv = a Af (a\v\) = a x M x 2 . (16) 

For a reheat temperature much smaller than M$, must be small. 
From Eq. (j|), the reheat temperature in terms of ax and Mx is T R h ~ 
a 1 ^ 2 M^Mpi. For = 10 13 GeV, a^ must be smaller than of order 
1CP 13 . On the other hand, ax may be as large as of order unity, or it may 
be small also. 

It is also convenient to work with dimensionless quantities that can 
absorb the effect of expansion of the universe. This may be accomplished 
with the definitions 



R 



X 



pxM^a 3 



(17) 



p R a , = 

It is also convenient to use the scale factor, rather than time, for the inde- 
pendent variable, so one may define a variable x = aM^. With this choice 
the system of equations can be written as (prime denotes d/dx) 
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(18) 



The constants c\, C2, and C3 are given by 
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3_Mpi 
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(19) 



Xeq is the equilibrium value of X , given in terms of the temperature T as 
(assuming a single degree of freedom for the X species) 
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EQ 
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(20) 



The temperature depends upon R and g*, the effective number of degrees 
of freedom in the radiation: 
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(21) 



It is straightforward to solve the system of equations in Eq. ( |18| ) with 
initial conditions at x — xi of R(xj) — X(xj) = and ^(xj) — $7. It 
is convenient to express p$(x — xj) in terms of the expansion rate at xi, 
which leads to 

3 Mp, H? o 

4> <t> 

The numerical value of xi is irrelevant. 

Before numerically solving the system of equations, it is useful to con- 
sider the early-time solution for R. Here, early times means H 3> T^, 
i.e., before a significant fraction of the comoving coherent energy density 
is converted to radiation. At early times $ ~ $7, and R ~ X ~ 0, so the 
equation for R' becomes R' — cix 3 ^ 2 ^ 1 / 2 . Thus, the early time solution 
for R is simple to obtain: 
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(23) 



Now express T in terms of R to yield the early-time solution for T: 
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(24) 



Thus, T has a maximum value of 
Tmax I 12 



0.77 



0.77a 




(25) 
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Figure 4. The evolution of energy densities and T/Mx as a function of 
the scale factor. Also shown is X/Xeq- 

which is obtained at x/xj = (&/?>) 2 / b = 1.48. It is also possible to express 
u§ in terms of Tjhj and obtain 

For an illustration, in the simplest model of chaotic inflation Hi ~ Ma, 
with Ma, ~ 10 13 GeV, which leads to T M ax/T rh ~ 10 3 (200/g*) 1 / 8 for 
T RH = 10 9 GeV. 

We can see from Eq. ( p3| ) that for i/a;/ > 1, in the early-time regime 
T scales as a -3 / 8 , which implies that entropy is created in the early-time 
regime [[H] . So if one is producing a massive particle during reheating it is 
necessary to take into account the fact that the maximum temperature is 
greater than T R n, and that during the early-time evolution, T oc a -3 / 8 . 

An example of a numerical evaluation of the complete system in Eq. 
( |l8|) is shown in Fig. || (from |p8| ). The model parameters chosen were 
Ma, = 10 13 GeV, a^ = 2x 1CT 13 , M x = 1.15 x 10 12 GeV, a x = 1(T 2 , and 
g* = 200. The expansion rate at the beginning of the coherent oscillation 
period was chosen to be Hi = M^. These parameters result in T R n = 
10 9 GeV and n x h 2 = 0.3. 

Figure ^ serves to illustrate several aspects of the problem. Just as ex- 
pected, the comoving energy density of <fi (i.e., a 3 pA,) remains roughly con- 
stant until ~ H, which for the chosen model parameters occurs around 
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a/aj ~ 5 x 10 8 . But of course, that does not mean that the tempera- 
ture is zero. Notice that the temperature peaks well before "reheating." 
The maximum temperature, Tmax = 10 12 GeV, is reached at a/ai slightly 
larger than unity (in fact at a/aj — 1.48 as expected), while the reheat 
temperature, Tjih = 10 9 GeV, occurs much later, around a/ai ~ 10 8 . Note 
that Tmax — W 3 Trh in agreement with Eq. (|26|). 

From the figure it is clear that X <C Xeq at the epoch of freeze out 
of the comoving X number density, which occurs around ajaj ~ 10 2 . The 
rapid rise of the ratio after freeze out is simply a reflection of the fact that 
X is constant while Xeq decreases exponentially. 

A close examination of the behavior of T shows that after the sharp ini- 
tial rise of the temperature, the temperature decreases as a -3 / 8 [as follows 
from Eq. (0)] until H ~ T^, and thereafter T oc a -1 as expected for the 
radiation-dominated era. 

For the choices of M^, a,p, g*, and ax used for the model illustrated in 
Fig. |, Q x h 2 = 0.3 for M x = 1-15 x 10 12 GeV, in excellent agreement with 
the mass predicted by using an analytic estimate for the result |p8| 

^=«(fr ,2 (^)'- ™ 

Here again, the results have also important implications for the con- 
jecture that ultra-high cosmic rays, above the Greisen-Zatsepin-Kuzmin 
cut-off of the cosmic ray spectrum, may be produced in decays of super- 
heavy long-living particles [|(], [32]] . In order to produce cosmic 
rays of energies larger than about 10 13 GeV, the mass of the A-particles 
must be very large, Mx ~ 10 13 GeV and their lifetime tx cannot be much 
smaller than the age of the Universe, tx ~ 10 10 yr. With the smallest 
value of the lifetime, the observed flux of ultra-high energy cosmic rays will 
be reproduced with a rather low density of A-particles, Qx ~ 10~ 12 . It 
has been suggested that X-particles can be produced in the right amount 
by usual collisions and decay processes taking place during the reheating 
stage after inflation if the reheat temperature never exceeded Mx f33]- 
Again, assuming naively that that the maximum number density of a mas- 
sive particle species X produced after inflation is suppressed by a factor 
of (Mx /Trh) 3 ^ 2 exp(— Mx/Trh) with respect to the photon number den- 
sity, one concludes that the reheat temperature Trh should be in the range 
10 11 to 10 15 GeV [^0|. This is a rather high value and leads to the gravitino 
problem in generic supersymmetric models. This is one reason alternative 
production mechanisms of these superheavy X-particles have been pro- 
posed p3L p3, p4[. However, our analysis show that the situation is much 



more promising. Making use of Eq. (27), the right amount of A-particles 
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to explain the observed ultra-high energy cosmic rays is produced for 




where it has been assumed that (<r\v\) ~ Mj{ 2 - Therefore, particles as 
massive as 10 15 GeV may be generated during the reheating stage in abun- 
dances large enough to explain the ultra-high energy cosmic rays even if 
the reheat temperature satisfies the gravitino bound. 

4-3. Production During Preheating 

Another way to produce wimpzillas after inflation is in a preliminary stage 
of reheating called "preheating" |2^] , where nonlinear quantum effects may 
lead to an extremely effective dissipational dynamics and explosive particle 
production. 

Particles can be created in a broad parametric resonance with a fraction 
of the energy stored in the form of coherent inflaton oscillations at the 
end of inflation released after only a dozen oscillation periods. A crucial 
observation for our discussion is that particles with mass up to 10 15 GeV 
may be created during preheating Q , and that their distribution 

is nonthermal. If these particles are stable, they may be good candidates 
for WIMPZILLAS [|7| . 

The main ingredient of the preheating scenario introduced in the early 
1990s is the nonperturbative resonant transfer of energy to particles in- 
duced by the coherently oscillating inflaton fields. It was realized that 
this nonperturbative mechanism can be much more efficient than the usual 
perturbative mechanism for certain parameter ranges of the theory p2| . 
The basic picture can be seen as follows. Suppose there is a scalar field X 
with a coupling g 2 <p 2 X 2 where 4> is a homogeneous classical inflaton field. 
The mode equation for X field then can be written in terms of a redefined 
variable \k = X^o?/ 2 as 

Xk (t) + [A + 2q cos(2t)} Xk (t) = (29) 

where A depends on the energy of the particle and q depends on the inflaton 
field oscillation amplitude. When A and q are constants, this equation is 
usually referred to as the Mathieu equation which exhibits resonant mode 
instability for certain values of A and q. In an expanding universe, A and 
q will vary in time, but if they vary slowly compared to the frequency of 
oscillations, the effects of resonance will remain. If the mode occupation 
number for the X particles is large, the number density per mode of the 
X particles will be proportional to |xfc| 2 - If A and q have the appropriate 
values for resonance, \k will grow exponentially in time, and hence the 
number density will attain an exponential enhancement above the usual 

16 



perturbative decay. This period of enhanced rate of energy transfer has 
been called preheating primarily because the particles that are produced 
during this period have yet to achieve thermal equilibrium. 

This resonant amplification leads to an efficient transfer of energy from 
the inflaton to other particles which may have stronger coupling to other 
particles than the inflaton, thereby speeding up the reheating process and 
leading to a higher reheating temperature than in the usual scenario. An- 
other interesting feature is that particles of mass larger than the inflaton 
mass can be produced through this coherent resonant effect. This has been 
exploited to construct a baryogenesis scenario p5j in which the baryon 
number violating bosons with masses larger than the inflaton mass are cre- 
ated through the resonance mechanism. A natural variation on this idea is 
to produce WIMPZILLAS by the same resonance mechanism. 

Interestingly enough, what was found |}7| is that in the context of a 
slow- roll inflation with the potential V(<f>) = m 2 ^(j> 2 /2 with the inflaton 
coupling of g 2 (f) 2 X 2 /2, the resonance phenomenon is mostly irrelevant to 
wimpzilla production because too many particles would be produced if 
the resonance is effective. For the tiny amount of energy conversion needed 
for wimpzilla production, the coupling g 2 must be small enough (for a 
fixed Mx) such that the motion of the inflaton field at the transition out 
of the inflationary phase generates just enough nonadiabaticity in the mode 
frequency to produce wimpzillas . The rest of the oscillations, damped 
by the expansion of the universe, will not contribute significantly to wim- 
pzilla production as in the resonant case. In other words, the quasi- 
periodicity necessary for a true resonance phenomenon is not present in 
the case when only an extremely tiny fraction of the energy density is con- 
verted into wimpzillas. Of course, if the energy scales are lowered such 
that a fair fraction of the energy density can be converted to WIMPZILLAS 
without overdosing the universe, this argument may not apply. 

The main finding of a detailed treatment J37j] is that WIMPZILLAS with 
a mass as large as 10 3 _ff/, where Hj is the value of the Hubble expansion 
rate at the end of inflation, can be produced in sufficient abundance to be 
cosmologically significant today. 

If the wimpzilla is coupled to the inflaton ^ by a term g 2 (j) 2 X 2 /2, then 
the mode equation in Eq. [l^ is now changed to 

cj 2 + k 2 + (M 2 x +g 2 2 )a 2 , (30) 

again taking £ = 1/6. 

The procedure to calculate the wimpzilla density is the same as in 



Section 4.1. Now, in addition to the parameter Mx/Hi, there is another 
parameter gMpi/Hj. Now in large-field models Hi ~ 10 13 GeV, so Mpi/Hi 
might be as large as 10 6 . The choice of g = 10~ 3 would yield gMpi/Hi = 
10 3 . 
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Figure 5. A graph of n x h 2 versus M x /Hi for gMpi/Hj = 10 6 . The solid 
curve is a numerical result, while the dashed and dotted curves are analytic 
approximations. 




Figure 6. An illustration of the nonmonotonic behavior of the particle 
density produced with the variation of the coupling constant. The value of 
Mx /Hj is set to unity. 
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Fig. H (from |j7|) shows the dependence of the wimpzilla density upon 
Mx/Hj for the particular choice gMpi/Hj — 10 6 . This would correspond 
to g ~ 1 in large-field inflation models where Mpi/Hi = 10 6 , about the 
largest possible value. Note that Qx ~ 1 obtains for Mj/ff/ « 10 3 . The 
dashed and dotted curves are two analytic approximations discussed in |37j , 
while the solid curve is the numerical result. The approximations are in 
very good agreement with the numerical results. 

Fig. ^ (also from |37| ) shows the dependence of the WIMPZILLA density 
upon gMpi/Hj. For this graph Mx/Hi was chosen to be unity. This 
figure illustrates the fact that the dependence of Uxh 2 on gMpi/Hj is not 
monotonic. For details, see p7|. 

4-4- Production in Bubble Collisions 

WIMPZILLAS may also be produced in theories where inflation is com- 
pleted by a first-order phase transition |38[ , in which the universe exits 
from a false- vacuum state by bubble nucleation |59| . When bubbles of true 
vacuum form, the energy of the false vacuum is entirely transformed into 
potential energy in the bubble walls. As the bubbles expand, more and 
more of their energy becomes kinetic as the walls become highly relativis- 
tic. 

In bubble collisions the walls oscillate through each other Gel and their 
kinetic energy is dispersed into low-energy scalar waves po|,^l|. We are 
interested in the potential energy of the walls, Mp = Airr/R 2 , where rj is the 
energy per unit area of a bubble wall of radius R. The bubble walls can be 
visualized as a coherent state of inflaton particles, so the typical energy E 
of the products of their decays is simply the inverse thickness of the wall, 
E ~ A -1 . If the bubble walls are highly relativistic when they collide, 
there is the possibility of quantum production of nonthermal particles with 
mass well above the mass of the inflaton field, up to energy A^ 1 = "fM^, 
with 7 the relativistic Lorentz factor. 

Suppose for illustration that the wimpzilla is a fermion coupled to the 
inflaton field by a Yukawa coupling g<fiXX. One can treat cj> (the bubbles 
or walls) as a classical, external field and the wimpzilla as a quantum 
field in the presence of this source. The number of WIMPZILLAS created in 
the collisions from the wall potential energy is Nx ~ fxMp/Mx, where 
fx parametrizes the fraction of the primary decay products in WIMPZIL- 
LAS. The fraction fx will depend in general on the masses and the cou- 
plings of a particular theory in question. For the Yukawa coupling g, it is 



fx — g ln(-yM^/2M x ) H |42 ■ wimpzillas may be produced in bubble 
collisions out of equilibrium and never attain chemical equilibrium. Even 
with TpH as low as 100 GeV, the present wimpzilla abundance would be 
fix ~ 1 if g ~ lO^a 1 / 2 . Here a -1 <§; 1 is the fraction of the bubble energy 
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at nucleation in the form of potential energy at the time of collision. This 
simple analysis indicates that the correct magnitude for the abundance 
of wimpzillas may be naturally obtained in the process of reheating in 
theories where inflation is terminated by bubble nucleation. 



5. Conclusions 

In this talk I have pointed out several ways to generate nonthermal dark 
matter. All of the methods can result in dark matter much more massive 
than the feeble little weak-scale mass thermal relics. The nonthermal dark 
matter may be as massive as the GUT scale, truly in the WIMPZILLA range. 

The mass scale of the wimpzillas is determined by the mass scale 
of inflation, more exactly, the expansion rate of the universe at the end 
of inflation. For large-field inflation models, that mass scale is of order 
10 13 GeV. For small-field inflation models, it may be less, perhaps much 
less. 

The mass scale of inflation may one day be measured! In addition to 
scalar density perturbations, tensor perturbations are produced in infla- 
tion. The tensor perturbations are directly proportional to the expansion 
rate during inflation, so determination of a tensor contribution to cosmic 
background radiation temperature fluctuations would give the value of the 
expansion rate of the universe during inflation and set the scale for the 
mass of the wimpzilla. 

Undoubtedly, other methods for wimpzilla production will be devel- 
oped. But perhaps even with the present scenarios one should start to 
investigate methods for wimpzilla detection. While wimpy wimps must 
be color singlets and electrically neutral, wimpzillas may be endowed with 
color and electric charge. This should open new avenues for detection and 
exclusion of wimpzillas. 

The lesson of the talk is illustrated in Fig. |?|. wimpzillas may surprise 
and be the dark matter, and we may learn that size does matter! 
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Figure 7. Dark matter may be much more massive than usually assumed, 
much more massive than wimpy wimps, perhaps in the wimpzilla class. 



References 

[1] B. Fuchs, "DARK98 Proceedings of the Second International Conference 
on Dark Matter in Astro and Particle Physics," eds. H V Klapdor- 
Kleingrothaus and L. Baudis, (Institute of Physics Publishing, Bristol 
and Philadelphia, 1999). 

[2] G. Evrard, "DARK98 Proceedings." 

[3] See, e.g., A. Dekel, in Proceedings of the 3rd ESO-VLF Workshop on 
"Galaxy Scaling Relations: Origins, Evolution and Applications," ed. L. 
da Costa (Springer). 



[41 See, e.g., J. A . Fyson, G. P. Kochanski, and I. P.. Dell'Antonio, astro 



ph/9801195 



[5] S. Sarkar, "DARK98 Proceedings." 
[6] G. Boerner, "DARK98 Proceedings." 

[7] E. W. Kolb and M. S. Furner, The Early Universe, (Addison- Wesley, Menlo 
Park, Ca., 1990). 

[8] K. Griest and M. Kamionkowski, Phys. Rev. Lett. 64, 615 (1990). 

[9] B. Cabrera, "DARK98 Proceedings." 
[10] I. Liubarsky, "DARK98 Proceedings." 
[11] R. Bernabei, "DARK98 Proceedings." 
[12] Y. Ramachers, "DARK98 Proceedings." 
[13] L. Baudis, "DARK98 Proceedings." 

21 



[14] 
[15] 
[16] 
[17] 
[18] 
[19] 
[20] 
[21] 
[22] 



J. Ellis, J.L. Lopez and D.V. Nanopoulos, Phys. Lett. B247, 257 (1990). 



See, for instance, G. F. Giudice and R. Rattazzi, hep-ph/9801271 
S. Raby, Phys. Rev. D 56, (1997). 



K. Hamaguchi, Y. Nomura and T. Yanagida, hep-ph/9805346 



K. Benakli, J. Ellis and D.V. Nanopoulos, |hep-ph/980333; . 
D. V. Nanopoulos, "DARK98 Proceedings." 

V. A. Kuzmin and V. A. Rubakov, Phys. Atom. Nucl. 61, 1028 (1998). 
tiep-ph/9804285| . 



M. Birkel and S. Sarkar, 



A. Kofman, A. D. Linde and A. A. Starobinsky, Phys. Rev. Lett. 73, 
3195 (1994); S. Yu. Khlebnikov and I. I. Tkachev, Phys. Rev. Lett. 77, 
219 (1996); Phys. Lett. B390, 80 (1997); Phys. Rev. Lett. 79, 1607 
(1997); Phys. Rev. D 56, 653 (1997); G. W. Anderson, A. Linde and 
A. Riotto, Phys. Rev. Lett. 77, 3716 (1996); see L. K ofman, The origin 



of matter in the Universe: reheating after inflation, astro-ph/9605155 



UH-IFA-96-28 preprint, 16pp., to appear in Relativistic Astrophysics: A 
Conference in Honor of Igor Novikov's 60th Birthday, eds. B. Jones and 
D. Markovic for a more recent review and a collection of references; see 
also L. Kofman, A. D. Linde and A. A. Starobinsky, Phys. Rev. D 56, 
3258 (1997); J. Traschen and R. Brandenberger, Phys. Rev. D 42, 2491 
(1990); Y. Shtanov, J. Traschen, and R. Brandenberger, Phys. Rev. D 
51, 5438 (1995). 



[23] 
[24] 
[25] 

[26] 

[27] 
[28] 
[29] 
[30] 



[31] 
[32] 

[33] 
[34] 



D. J. H. Chung, Edward W. Kolb, and A. Riotto, [hep-ph/980223£ . 
S. Fulling, Gen. Rel. and Grav. 10, 807 (1979). 

N. D. Birrell and P. C. W. Davies, Quantum Fields in Curved Space (Cam- 
bridge University Press, Cambridge, 1982). 

D. M. Chitre and J. B. Hartle, Phys. Rev. D 16, 251 (1977); D. J. Raine 
and C. P. Winlove, Phys. Rev. D 12, 946 (1975); G. Schaefer and H. 
Dehnen, Astron. Astrophys. 54, 823 (1977). 

N. D. Birrell and P. C. W. Davies, J. Phys. A:Math. Gen. 13, 2109 (1980). 



D. J. H. Chung, E. W. Kolb, and A. Riotto, hep-ph/9809453 



D. H. Lyth and A. Riotto, tiep-ph/9807278 



J. Ellis, J. Kim and D. V. Nanopoulos, Phys. Lett. B145, 181 (1984); L. 
M. Krauss, Nucl. Phys. B227, 556 (1983); M. Yu. Khlopov and A. D. 
Linde, Phys. Lett. 138B, 265 (1984); J. Ellis, D. V. Nanopoulos, and S. 
Sarkar, Nucl. Phys. B461, 597 (1996). 

R. J. Scherrer and M. S. Turner, Phys. Rev. D31, 681 (1985). 

V. Berezinsky, M. Kachelriess and A. Vilenkin, Phys. Rev. Lett. 79, 4302 
(1997). 



V. Kuzmin and I. I. Tkachev, tiep-ph/9802304 



D. J. H. Chung, E. W. Kolb and A. Riotto, hep-ph/9805473 



22 



[35] E. W. Kolb, A. D. Linde and A. Riotto, Phys. Rev. Lett. 77, 4290 (1996). 

[36] E. W. Kolb, A. Riotto and I. I. Tkachev, Phys. Lett. B423, 348 (1998). 

[37] D. J. H. Chung, [hep-ph/980948£ . 

[38] D. La and P. J. Steinhardt, Phys. Rev. Lett. 62, 376 (1989). 

[39] A. H. Guth, Phys. Rev. D 23, 347 (1981). 

[40] S. W. Hawking, I. G. Moss and J. M. Stewart, Phys. Rev. D 26, 2681 (1982). 

[41] R. Watkins and L. Widrow, Nucl. Phys. B374, 446 (1992). 

[42] A. Masiero and A. Riotto, Phys. Lett. B289, 73 (1992). 



23 



